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O ! Abstract 

We explicate a procedure to solve general linear differential equations, which con- 
nects the desired solutions to monomials x x of an appropriate degree A. In the process 
the underlying symmetry of the equations under study, as well as that of the solutions 
are made transparent. We demonstrate the efficacy of the method by showing the 
common structure of the solution space of a wide variety of differential equations viz. 
Hermite, Laguerre, Jocobi, Bessel and hypergeometric etc. We also illustrate the use 
of the procedure to develop approximate solutions, as well as in finding solutions of 
(Jh many particle interacting systems, 

-i— > ' 

1 Introduction 

As is evident, the basic building blocks of polynomials and functions are the monomi- 
als. Superposition of finite and infinite number of these constituents, with appropriate 
coefficients, respectively leads to various polynomials and functions. Well known clas- 
sical orthogonal polynomials and functions like Bessel, Airy etc., manifest as solutions 
of Schrodinger equation in quantum mechanics, in the study of waves, in optics and 
' electromagnetism. 
{Sj ■ Apart from finding solutions of various differential equations, considerable effort 

has been put into understanding the symmetry properties of these equations, which 
throws light on the structure of the solution space [1]. Here, we illustrate a procedure 
of solving linear differential equations, which leads to an exact mapping between the 
$_i , desired polynomials and functions and the space of monomials [2, 3]. In this method, 

the Euler operator (EO) D = x-^ plays a crucial role. It makes essential use of the 
two properties of the EO. 

1. The monomials are the eigenfunctions of the EO: 

x-^-x x = Dx x = Xx x . (1) 
ax 

2. The Euler operator 'measures' the degree of an operator, irrespective of its 
constituents, 

[D, O d ] = dO d . (2) 
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In the above O d is an operator of 'degree' d e.g., x 2 and x 3 -^ have degree 2, while 
^7 and J^- + ^ have degree -2. 

In the following, we start with the familiar Hermite, Laguerre and confluent hyper- 
geometric equations and make use of their symmetry properties to obtain solutions in a 
form, which exhibits the mapping between the respective polynomials and the space of 
monomials. We then explicate a more general method, applicable to linear differential 
equations of arbitrary order and obtain the desired solutions making use of the above 
mentioned properties of the Euler operator. The procedure is explicated through so- 
lution of the oscillator problem. Subsequently, earlier solved examples are connected 
with the new approach, considering Hermite differential equation as the example. After 
illustrating the solution space of the hyper geometric equation, we tabulate the solu- 
tions and the structure of their respective differential equations for all the functions 
and polynomials, commonly encountered in physics literature. The case of the so called 
quasi exactly solvable Schrodinger equation is then studied, for which only a part of 
the spectrum is analytically obtainable. This example shows the use of the present 
approach for finding the exact solution, as well as in the development of approximation 
schemes. We then proceed to many-body interacting Calgero-Sutherland systems and 
show how our procedure can be profitably explored for finding out the eigenspectra 
and eigenfunctions. We then conclude with directions for future work. 



2 Connecting Polynomials with Monomials 

We start with the well studied Hermite differential equation: 

[~-^ + (D-n)]H n (x) = 6H n (x) = 0, (3) 

which has seen the use of traditional series solution method to the algebraic approach of 
raising and lowering operators, for finding the solution [4]. It is worthwhile to observe 
that the above equation has the Euler operator and a degree d=-2 operator — = P, 
apart from a harmless constant. The fact that [D,P] = —IP can be profitably used 
through the Baker-Campbell-Hausdorff (BCH) formula: 

e A Be- A = B + [A,B)+ l - [A, [A, B}} + (4) 

in order to map the differential operator O onto (D-n). 
For this purpose, it can be straightforwardly checked that 

id 2 Id 2 _i d 2 Id 2 Id 2 

e ^ [ -2d^ 2 + D ~ n]e ~^ = -2dx^ + D - n+ 4 [ ^ D] = (D - n) - (5) 

This immidiately suggests to factorize the Hermite polynomial in the form 

1 d 2 _i d 2 _i d 2 

[--^-2 + D - n]e 4dx?ip(x)=e 4 dx 2 ~[D-n](p(x) = 0, (6) 

which implies, 

[D - n]<p(x) = 0. (7) 
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The fact that the monomials are eigenfunctions of the Euler operator D leads to <f>(x) = 

c n x n and hence the Hermite polynomial becomes, H n (x) = c n e 4 a^x™ [5]. The above 
is the exact connection between the monomial x n and the polynomial H n (x); the value 
of c n is to be determined from normalization condition. The discerning reader can 
visualize that the structure of Laguerre differential equation, 

^-"-< a+1 >s-^ i »w = °- (8) 

and confluent hyper geometric equation, 

[x^+a-x^-7^]*(a >7 ,x) = > (9) 

are similar to the Hermite case. A similar mapping connects the above differ- 
ential operators to (D — n) and hence the solutions can be written in the form, 
L a n = c n exp[-x-£^ - (a + l)£]x n and $(a, 7, x) = c^ a exp[-x-£? - l-^]x~ a . 

This procedure can be used to identify the generating function and the algebraic 
structure of the solution space [6] . We refer the interested reader to Ref [2] for more 
details, as also to Ref [7] for construction of coherent states for the dynamical systems, 
associated with the harmonic oscillator and the Coulomb problem. 

It is evident that, if the differential operator O contains an operator of definite 
degree, apart from the Euler operator and the constant term, the aforementioned ap- 
proach can connect the solution space to monomials through an exponential mapping. 
It may not be the case in general. Keeping this in mind, in the next section we develop 
a more general approach connecting the solution space with monomials for general 
equations. 



3 A simple approach to familiar differential 
equations 

For simplicity, we will first consider the case of single variable linear differential equa- 
tions (LDEs) and point out its multi-variate generalization later. A single variable 
LDE, as will become clear from the examples of later sections, can be cast in the form 

[F(D) + P(x,d/dx)]y(x) = 0, (10) 

where, F(D) = J2n=^oo a nD n and a„'s are some parameters; P(x,d/dx) can be an 
arbitrary polynomial function of x, d/dx and other operators. The solution to Eq. (10) 
can be written as [3], 

{00 
£(-!) 
m=0 

provided, F(D)x x = 0; here C\ is constant. It is easy to see that, the operator 
1/F(D) is well defined in the above expression and will not lead to any singularity, 
if P(x,d/dx) does not contain any degree zero operator. We note that D itself is an 
operator of degree zero. 



F(D) 



P(x,d/dx) 



(11) 
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The proof of Eq.(ll) is straightforward and follows by direct substitution [3]. Al- 
ternatively, since F(D)x x = 0, equating F(D)x x modulo C\ and Eq. (10), one finds, 



[F(D) + P(x, d/dx)} y(x) = C\F(D)x x 
Rearranging the above equation in the form 



F(D) 



1 + 



1 



F(D) 

and cancelling F(D), we obtain 

1 



P(x, d/dx) 



1 + 



This yields 



F(D) 

y(x) = c x 



P(x,d/dx) 



y(x) = C x F(D)x ) 



y(x) = C\x > 



1 



1 + F{D] P ( X 'd/dx) 

which can be cast in the desired series form: 

( _ ir _1 

m=0 



(12) 
(13) 

(14) 
(15) 



y(x) = C x J2(-l) 



F{D) 



P(x, d/dx) 



It is explicit that, the above procedure connects the solution y(x) to the space of the 
monomials x x . The generalization of this method to a wide class of many- variable 
problems is immediate. Using the fact that, F(D)X X = has solutions, in the space 
of monomial symmetric functions [8], where D = = J2i x i7f^i the solutions of 

those multi-variate DEs, which can be separated into the form given in Eq. (10), can 
be solved like the single variable case. As will be seen later, this procedure enables one 
to solve a number of correlated many-body problems. 

For illustration, we consider the harmonic oscillator problem. The Schrodinger 
eigenvalue equation (in the units, H=u=m=l) 



dx 2 



+ (2E n - x 2 ) 



Ipr. 







can be written in the form given in Eq. (10), after multiplying it by x 2 : 

(D- l)D + x 2 (2E n -x 2 )] 4> n = . 



(16) 



(17) 



Here, F[D) = (D - 1)D and the condition F(D)x x = yields, A = or 1. Using Eq. 
(11), the solution for A = is, 



V>o = Co EH 

lm=0 



■(x 2 (2E -x 2 )) 



Co 



1 



(D-1)D 

[2E ] 2 , (2\ + [2E ] 2 ) 4 (4! + (2!) 2 [2S ] + (2!)[2i? ] 3 ; 



2! 



-x + 



4! 



-x 



2!6! 



x 6 + • • ( L8.) 



Here ipo is an expansion in powers of x, whose coefficients are polynomials in Eq. The 
above series can be written in a closed, square integrable form, Co exp(— x 2 /2), only 
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when Eq = 1/2. Analogously, A = 1, yields the first excited state. To find the n th 
excited state, one has to differentiate the Schrodinger equation (n — 2) number of times 
and subsequently multiply it by x n to produce a F(D) = x n £^t = ]Xi=q{D — I), and 
proceed in a manner similar to the ground state case. 

It is clear that, our procedure yields a series solution, where additional conditions 
like square integrability has to be imposed to obtain physical eigenfunctions and their 
corresponding eigenvalues. Once, the ground state has been identified and for those 
cases, where ip(x) = ipoP(x), where P (x) is a polynomial, one can effortlessly obtain 
the polynomial part, as will be shown below. Proceeding with the harmonic oscillator 
case and writing 

x 2 

ip a (x) = exp( — —)H a (x) , 



one can easily show that H a satisfies 



D-a- 



2dx 2 



H a (x) = 



where a = E n — 1/2. The solution of the DE 

oo 

H a {x) = c a Y,{-vr 



m=0 



1 d 2 



(D - a) 2 dx 2 



(19) 



(20) 



yields a polynomial only when a is an integer, since the operator d?/dx 2 reduces 
the degree of x a by two, in each step. Setting a = n in Eq. (19), we obtain the 
Hermite DE and E n = (n + 1/2) as the energy eigenvalue. Below, we give the algebraic 
manipulations required to cast the series solution of Eq. (20) into a form, not very 
familiar in the literature. For the Hermite DE, F(D) = D—n and P(x, d/dx) - 
the condition F(D)x x = yields A = n, hence, 



l d 2 
2E 1 ' 



m=0 



(D 



n) 2 dx 2 



Using, [D , (d 2 /dx 2 )} = —2(d 2 /dx 2 ) and making use of the fact 



dse-< D ~^ 



(D 



n) 



(21) 



(22) 



we can write, 



1 1 



d 2 



2 (D - n) dx 2 



1 d 2 



1 



1 1 



d 2 



2{D-n) dx 2 
Hence in general, 



d 2 



2 (D - n) dx 2 



d 2 



2 dx 2 {D — n- 2) 

2 



2^2 



1 



IP 



n 



4) (D 



n 



2) 



(23) 



2 (D - n) dx 2 




(24) 
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Substituting Eq. (24) in Eq. (21), we obtain, 



m=0 \ / 



= C n exp 



4 cte 2 



(25) 



a result, not commonly found in the literature. The arbitrary constant C n is chosen to 
be 2™, so that the polynomials obtained can match with the standard definition [9]. 
Likewise, the solution to the hyper geometric DE 



2 d , n -,\ d n d 2 d 
x — ^ + (a + p + 1) x— + ap - x—^ - 7— 
dx z dx dx z dx 



F(a,f3; T ,x)=0 



(26) 



can be written as, 



F(a,(3; T ,x) = (-l) 



r( 7 -/3)r(«) 



exp 



(D 



a 



, dx 2 ^ ^ dx 



. x 



(27) 



For convenience a table has been provided at the end, which lists a number of commonly 
encountered DEs and the novel exponential forms of their solutions. It is worth point- 
ing out that, the solution of hypergeometric DEs is a polynomial solution, provided 
either a or f3 are negative integers. It should be noticed that, unlike the conventional 
expressions, the monomials in the above solutions are arranged in decreasing powers 
of x. 

Multiplying the hypergeometric DE with x yields two roots, A = 0, 1 — 7 and A = 
solution gives rise to the well-known, Gauss hypergeometric series. Since a number 
of quantum mechanical problems can be related to confluent and hypergeometric DEs 
[10], we hope that the novel expressions given above and in the table will find physical 
applications. 



4 Quasi-exactly solvable problems 

This section is devoted to the study of QES problems [11]. These problems are in- 
termediate to exactly and non-exactly solvable quantum potentials, in the sense that, 
only a part of the spectrum can be determined analytically. These potentials have 
attracted considerable attention in recent times, because of their connection to various 
physical problems [12]. 

We illustrate our procedure, through the sextic oscillator in the units H = 2m = 
lj = 1, for the purpose of comparison with the standard literature. The corresponding 
eigenvalue equation is given by: 

if,(x) = EiP(x) . (28) 

Asymptotic analysis suggests a trial wave function of the form, 

il>{x) = exp(-6x 4 )^(x) , (29) 



d 2 

2 6 
- -r^ + OiX +JX 
dx z 
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which leads to, 



dx 2 



2^x 3 — + (a + 3^)x 2 



ijj(x) = E^(x) , 



(30) 



where x 6 term has been removed by the condition 166 2 = 7. One notices that, the 
operator O = (a + 3^/t~)x 2 + 2^fyx 3 d/dx increases the degree of ip(x) by two, if ^(x) 
is a polynomial. Confining ourselves to polynomial solutions and assuming that the 
highest power of the monomial in ip(x) is n, one obtains, 



a 



= 2n + 3 



(31) 



after imposing the condition that O does not increase the degree of the polynomial. 
This is the well-known relationship between the coupling parameters of the quasi- 
exactly solvable sextic oscillator [12]. Taking n = 4 and 7=1 for simplicity, and after 
multiplying the above equation with x 2 : 



D(D - 1) + Ex 2 + 8x 4 - 2x l 



d_ 

dx 



i){x) = 



we get, 



>m=0 



_J_( i2£o + 8i1 _ 2i ^ 



Modulo Co, the above series can be expanded as, 

„2 /t 



"00 0*0 



1 - E - + (E 2 - 16)- + (64^ - E*)- + 



(32) 



(33) 



(34) 



The monomials having degree greater than four vanish provided, £^0 = 0, ±8. It can be 
explicitly checked that, for these values of Eq, Eq. (32) is satisfied. The eigenfunctions 
corresponding to these three values are given by, 



and 



V>_ 8 (x) = exp(-— )[l + 4x 2 + 2x 4 ] 

x 4 2 
ifo(x) = exp(-— )[1 - -x 4 ] , 

x 4 

ip +8 (x) = exp(-— )[1 -4x 2 + 2x 4 ] 



(35) 
(36) 
(37) 



This procedure generalizes to a wide class of QES problems [13]. 

Below, we demonstrate the method of finding approximate [13] eigenvalues and 
eigenfunctions for non-exactly solvable problems, using the well studied anharmonic 
oscillator as the example : 



d 2 

— r-Q + ax<2 + P xA 
dx z 



En 



ip n (x) = 



Proceeding as before, V'o(^) can be written as, 

1 

(_1)'" 

1 ; (71. 

,m=0 



Mx) = Co{ £(-1) 



(D-1)D 



(Eqx 2 



ax 



/3x ( 



.1 



(38) 



(39) 
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which can be expanded as, 



Y> (x) = 1 - ^x 2 + i (2a + E 2 ) x 4 - I (24/3 - (14a£ + £ 3 )) x 6 + • • • (40) 

Although a number of schemes can be devised for the purpose of approximation, we 
consider the simplest one of starting with a trial function tpo(x) = exp(— fix 2 — vx A ) 
and matching it with ipo(x)- Comparison of the first three terms yields, 



E 

2! 4! ^ 4! ' 

and = A_(^o^o 3 ) . ( 41 ) 

p 3! 30 6! K ' 

The resulting cubic equation in energy, Eq — Ega = 3/3/2, leads to one real root and 
two complex roots. Choosing the real root on physical grounds, one obtains, 

r, 2l/3 « A 

^0 = ^- + ^ , (42) 

r 1 1/3 

where A = 40.5/3 + (1640.25^ 2 - 108a 3 ) 1 / 2 . The value of E , obtained in the 

weak coupling regime, matches reasonably well with the earlier obtained results [14]. 
An approximate tpo can be obtained from Eq. (40). One can easily improve upon the 
above scheme by taking better trial wave functions. Similar analysis can be carried 
out for the excited states. The above expansion of the wave function may be better 
amenable for a numerical treatment. For example, an accurate numerically determined 
energy value can lead to a good approximate wave function. 



5 Many-body interacting systems 

In this section, we will be dealing with correlated many-body systems, particularly 
of the Calogero-Sutherland [15] and Sutherland type [16]. These models have found 
application in diverse branches of physics like fluid flow, random matrix theory, novel 
statistics, quantum Hall effect and others [2]. 

We start with the relatively difficult Sutherland model, where the particles are 
confined to a circle of circumference L. The two-body problem treated explicitly below, 
straightforwardly generalizes to particles. The Schrodinger equation is given by (in 
the units h = m = 1) 



i 2 d 2 



7T 



1=1 1 



L 2 sin 2 [7r(xi - x 2 )/L] 



E x 



MM) = o 



(43) 



Taking, zj = e 2mx j/ L and writing ip\({zi}) = Yli,i^j z i ^^{ z i ~ z jY J\{{zi}), the above 
equation becomes, 



£ D 2 + ^-^(D, - D 2 ) + Eq — E\ 



i=l 



J\({zi}) = 



(44) 
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where, A = z t £-, E\ = 2(±) 2 E X , E = 2(±) 2 E and E = (f ) 2 /? 2 , is the ground- 
state energy. Here, J\{{zi}) is the polynomial part, which in the multivariate case is 
the well known Jack polynomial [8]. Here, A is the degree of the symmetric function 
and {A} refers to different partitions of A. J2iD 2 is a diagonal operator in the space 
spanned by the monomial symmetric functions, m^x}, with eigenvalues X)?=i A 2 • A 
monomial symmetric function is a symmetrized combination of monomials of definite 
degree. For example for two particle case, there are two monomial symmetric functions 
having degree two. These are 777.2,0 = x\ + x\ and ra^i = x\x 2 . Readers are referred 
to Ref. [9] for more details about various symmetric functions and their properties. 
Rewriting Eq. (44) in the form, 



E(A 2 - A 2 ) + -D 2 ) + E Q + Z A 2 - E\ 



z\ ~ z 2 

one can immediately show that, 



ln=0 



J x ({ Zi }) = , (45) 



E,(A 2 -A 2 ) \ r Zl -z 2 



xmx({ Zl }) . (46) 



For the sake of convenience, we define 

1 



5 = 
and Z = (3 



E^-A 2 

Z\ + Z 2 



z\ - z 2 



(D 1 -D 2 ) + E + J2>H-E^ 



(47) 



The action of S on mx({zi}) yields singularities, unless one chooses the coefficient of 
mx in Z mx({zi}) to be zero; this condition yields the eigenvalue equation 

Ex = E + J2(^i +/3[3-2i]Ai) 



Using the above, one can write down the two particle Jack polynomial as, 

EM^f) ( ^f <D '- D2) -? (3 - 2i)A ') 



n=0 



xmx({xi}) . (48) 



Starting from 777,2,0 = z\ + z|, it is straightforward to check that 
Zm 2fi = 4/5(2;i + z 2 f = 4/3m 2 

^ m2 >° = -4) ^ m lo^ = - 2 ^ m lo 

S n m 2fi = -2{J3) n m\ fl for n > 1 . 

Substituting the above result in Eq. (46), apart from C2, one obtains 



■h = m 2 ,o+ £(-in-2)(/3r Wo 



\n=l 
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\n=0 



m 2 n + 



1 + 



m 



1,0 



(49) 



which is the desired result. The above approach can be easily generalized to the N- 
particle case. 

Another class of many-body problems, which can be solved by the present approach 
is the Calogero-Sutherland model (CSM) and its generalizations. Proceeding along the 
line, as for the Sutherland model, one finds the eigenvalues and eigenfunctions for the 
CSM. 

The Schrddinger equation for the CSM in the previous units, is given by, 



2 N 



2^&^ + 2^ * + 2 ^ (,, - xtf 



En 



*Pn({Xi}) = , 



(50) 



where the wave function is of the form [15], 

lP n (x) = ^Pn{{Xi}) = ZGP n ({ Xi }) . 



(51) 



Here Z = Hi^i ~ Xjf, G = expf-^x? } , 9 2 = P(P ~ 1) and P n {{x t }) is a 
polynomial. After removing the ground state, the polynomial P n {{xi}) satisfies, 



— 2 -pY1 l _ — 



Pn({Xi}) = 



(52) 



where Eq = IN + 1/3 N(N — 1). Defining 







Id 2 1 



2 y dx i 

one can easily see, following the procedure adopted for the Hermite DE that, 

P n ({xi}) = C n e- A Wm {n} ({ Xi }) , 



(53) 



where m^ n }({xi}) is a monomial symmetric function of degree n. The corresponding 
energy is given by 

E n = E + n , (54) 
One can use the above procedure to solve many other interacting systems. 



6 Conclusions 

In conclusion, we have presented a novel scheme to treat exactly, quasi-exactly and 
non-exactly solvable problems, which also extends to a wide class of many-body inter- 
acting systems. The procedure can be used for the construction of the ladder operators 
for various orthogonal polynomials and the quantum systems associated with them. 
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The approximation scheme presented needs further refinement. It should be analyzed 
in conjunction with computational tools for finding its efficacy as compared to other 
methods. The many-body problems presented here have deep connection with diverse 
branches of physics and mathematics. The fact that the procedure employed for solving 
them, connects the solution space of the problem under study to the space of mono- 
mials, will make it useful for constructing ladder operators for the many-variable case. 
This will throw light on the structure of the Hilbert space of these correlated systems. 
Some of these questions are currently under study and will be reported elsewhere. 
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TABLE I. Some frequently encountered DEs and their novel solutions. 

Below, the differential equations from top to bottom, respectively, are Hermite, Laguerre, 
Legendre, Gegenbauer, Chebyshev Type I, Chebyshev Type II, Bessel, confluent hypergeometric 

and hypergeometric. 



Differential Equation 



F(D), D = x{d/dx) 



Solution 



*±-n-(a + l)±-x^ 



L« = 



x^ + (2\ + l)x±-n(2\ + n)-£ 
xChx) = 



x 



2 J$_ _j_ rp_d_ ^2 d 2 

X dx? X dx n dx'- 



dx* ' " ax ' ' ax* ~ 

-n(n + 2)-£ 



T n (x 
d 2 



} 2 £+ x £-» 2 + x ' 



= 

U n (x] 

= 



= 



(X 



>(a,7,x) 



= 



x(l-x)^ + ( 7 - [a + 



x 



(D — n) 

(D-n) 

(L> + ra + 1)(£> -n) 
(D + n + 2X)(D -n) 

{D + n){D-n) 
(D + n + 2)(D-n) 

{D + u){D-u) 

(D + a) 

(D + a)(D + f3) 



H n (x) = C„exp 



.1 d? 

idx 1 



.X" 



L%(x) = C n exp 



P n (x) = C„exp 
C*(x) = C„exp 



T n (x) = C„exp 
U n (x) = C n exp 



J±v(x) = C±^exp 
$(a,7,x) = C_Q,exp 



(a + 1) 



dx 



.X 



1 d 2 " 

2(D+n+l) IS 1 

1 d 2 ~ 
2(D+rt+2A) ch? 



.X" 



1 d 2 



2(D+rt) dx 1 
1 d 2 
2(D+n+2) di 2 " 



1 X 2 



.X 



.X 



2(D±i/) 



.X 



TV 



x dx^ I dx 



xF(a, (3, 7, x) = 



x exp 



F(a,/3, 7 ,x) = C_ (aj/3) x 



-G9,a) 



The solution to the DE [F(D) + P(x, d/dx)} y(x) = 0, is, 
V(x) = C A {E™=o(-l) m [ F ( 7JJ P(x,d/dx)] m }x A , provided F(D)x A = 0. 
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